Role of kinetic inductance in transport properties of shunted superconducting nanowires 
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Recently transport measurements have been carried out in resistively shunted long superconducting nanowires 
[M. W. Brenner et. ah, Phys. Rev. B 85, 224507 (2012)]. The measured voltage-current iV-I) characteristics 
was explained by the appearance of the phase slip centers in the shunted wire, and the whole wire was modeled 
as a single Josephson junction. The kinetic inductance of the long nanowires used in experiments is generally 
large. Here we argue that the shunted superconducting nanowire acts as a Josephson junction in series with an 
inductor The inductance depends on the length and the cross section of the wire. The inclusion of inductance 
in our analysis modifies the V-I curves, and increases the rate of switching from the superconducting state to 
the resistive state. The quantitative differences can be quite large in some practical parameter sets, and might be 
important to properly understand the experimental results. Our proposed model can be verified experimentally 
by studying the shunted superconducting nanowires of different lengths and cross sections. 

PACS numbers: 74.78.Na, 74.25.F-, 74.25.Sv, 74.40.-n 



I. INTRODUCTION 

Superconducting nanowires (SCNWs) of radius of several 
nanometers and length of hundreds of nanometers are under 
an intense current research for studying various interesting 
physical phenomena and potential technological applications 
including the fluctuations in superconductors, the supercon- 
ducting qubits and the current standards ' . The voltage-current 
{V-I) characteristics of the SCNWs are measured by several 
groups, and a hysteresis has been observed. Upon increas- 
ing the bias current, the SCNW switches to a resistive state 
at a critical current /jw. When one reduces the bias current 
from the resistive state, the nanowire remains in the resistive 
state until the current is below a retrapping current 1^. The 
V-I curves of a SCNW are very similar to those of an under- 
damped Josephson junction (JJ). 

By numerically solving the generalized time-dependent 
Ginzburg-Landau equations, it was shown in Refs. 2 and 3 
that the hysteresis in the SCNWs can be explained by allow- 
ing the occurrence of a phase slip center (PSC). The super- 
conducting phase in a PSC evolves with time according to 
the ac Josephson relation, and the superconductivity is sup- 
pressed, while the phase at the two sides of the PSC remains 
constanPI^. The PSC provides a mechanism for survival of 
superconductivity even in the presence of an electric field in- 
side the superconductor The SCNW with a PSC can be mod- 
eled as a JJ. When the SCNW is driven by a dc bias current, a 
hysteresis shows up in the observed V-I curves. However for 
a voltage biased case, the reported V-I curves are S-shapecP. 
For a review on the PSC, please check Ref. 2i 

The transport measurements are generally performed with 
the free-standing nanowires'^''^. The exchange of heat from 
the free-standing SCNWs is weak because the cooling occurs 
only at the two ends of the nanowire where they are attached 
to superconducting electrodes. The self-heating thus becomes 
important for these nanowires. The hysteretic V-I curves were 
obtained in a theoretical analysis by Tinkham et. a/.^by us- 
ing the measured temperature T dependent linear resistance 



R(T) as an input in solving the heat flow equation through 
the SCNWs. At some current interval (between I^ and 1^,^), 
there are two possible states with one being the heating dom- 
inated resistive state and the other superconductivity domi- 
nated state. A large fraction of the nanowire becomes nor- 
mal in the heating dominated resistive state while only a small 
fraction of the wire is normal in the superconductivity domi- 
nated state. Therefore the nanowire shows the hysteretic V-I 
curves. The self-heating model was further advanced in Ref. 
9 

The switching current /jw becomes stochastic in the pres- 
ence of thermal or/and quantum fluctuations. A distribution 
of is obtained in experiments by repeating the measure- 
ments of /sw. It has been demonstrated for the free-standing 
unshunted nanowires that the width of the distribution of /sw 
increases with a decreasing temperature** This counter- 
intuitive result was argued from the self-heating model. The 
unshunted SCNW is abounded with thermal phase slips at rel- 
atively high temperatures. The phase slip fluctuations induce 
resistance which causes Joule heating in the bulk, and the tem- 
perature of SCNW increases above the electrode temperature. 
When the electrode temperature is reduced, the self-heating 
caused by the dissipative phase slips becomes even more im- 
portant as a result of the reduced heat capacity and heat con- 
ductivity. It was also observed that the retrapping at 1^ is not 
stochastic for the unshunted SCNW within the experimental 
resolution. 

The V-I characteristics of the nanowires with an exter- 
nal shunt resistor have been measured last year". The V-I 
curves of a shunted SCNW changes dramatically in compari- 
son to that of an unshunted SCNW. While the mean switching 
current depends weakly on the external shunt resistance, the 
mean retrapping current increases with a decreasing shunt re- 
sistance. The mean retrapping current merges with the mean 
switching current at a threshold shunt resistance, and the hys- 
teresis in V-I curves vanishes. The distribution width of the 
switching current reduces with a decreasing shunt resistance 
or temperature in the low temperature region. All these obser- 
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vations cannot be explained by the self -heating model, which 
is only valid for the unshunted nanowires. The self-heating is 
reduced for a shunted nanowire because when there is a phase 
slip in the nanowire, the resistance of the wire is finite which 
causes a redistribution of the applied bias current in the shunt 
circuit. The authors of Ref. fTTlhave modeled the SCNW with 
a PSC as a JJ. The derived V-I curves and the distributions 
of the switching current based on their model of a JJ with an 
external shunt show reasonable qualitative agreement with the 
experimental results. 

The role of the rest part of the nanowire excluding the PSC 
at a weak link was completely neglected in the RCSJ type 
modeling of Ref. \n\ After a PSC develops in the SCNW, 
the bias current is redistributed between the nanowire and 
the shunt resistor Then an electric field is induced in the 
SCNW resulting from the time-dependent supercurrent. In 
this work, we show that the rest part of the SCNW apart from 
the PSC acts as a kinetic inductor due to the kinetic energy of 
Cooper pairs, which introduces a retardation in the dynamics 
of the superconducting phase. The value of the inductance de- 
pends on the length and the cross section of the wire. The ki- 
netic inductance affects the retrapping current, and also trans- 
forms the distribution of the switching current. Furthermore, 
the kinetic inductance reduces the self-heating effect in the 
nanowire. Thus, we conclude that the kinetic inductance of 
the SCNW cannot be neglected. 

The remaining part of the paper is organized as follows. In 
Sec. II, we derive the quantum circuit equations for our model. 
In Sec. Ill, the V-I curves without noise at zero temperature 
are calculated. In Sec. IV, the effect of the kinetic inductance 
on the rate of the switching from the superconducting state 
is investigated. A short discussion on our results and some 
comparisons with other related recent studies are presented in 
Sec. V. 



II. QUANTUM CIRCUIT AND NONLOCAL DAMPING 

We model the SCNW with a PSC as a JJ connected in se- 
ries with an inductor The resistance of the JJ is due to the 
phase slips, and can be measured experimentally. The induc- 
tor describes the rest part of the superconducting region ex- 
cluding the PSC. The inductance can be obtained from the 
London equation which states that the time-varying super- 
current density Js induces an electric field E inside the su- 
perconductor E = 4nA^diJs/c^ where A the London penetra- 
tion depth. Therefore, the kinetic inductance of the SCNW is 
Lk = V/d,h = 4nA^l/(c^A), where I, = J,A, I is the length 
of the wire and A is the cross section area of the nanowire. 
When the bias current is smaller than /jw, the whole bias cur- 
rent flows through the SCNW in the absence of phase slips. 
In the presence of phase slips in the nanowire, the resistance 
of the wire becomes finite and the current redistributes in the 
shunt circuit. The kinetic inductance Lk introduces a time 
scale for the redistribution ti - Lk/(^s +^psc) where Rp5c,Rs 
are the resistance of the PSC and the external shunt resistor re- 
spectively. As the phase slip disappears, the nanowire recov- 
ers to the superconducting state. The current flows fully again 
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FIG. L (color online) (a) A schematic of a dc current biased super- 
conducting nanowire with an external shunt resistor. The nanowire 
with an external shunt can be represented by an equivalent circuit 
shown in (b). The nanowire is modeled by a resistively capacitively 
shunted Josephson junction (RCSJ) in series with a kinetic inductor. 



into the nanowire, and the time scale for this redistribution is 

TP 

T"2 = I^k/Rs- In the experiment , Rpsc > so T2 > tj. Thus 
there is a sufficient spare time for the nanowire to cool down 
before the bias current redistributes into the wire. Therefore, 
the inductance reduces the heating effect in the SCNWs with 
a small external shunt resistor. The reduced self-heating ef- 
fect due to the kinetic inductance of the SCNWs was demon- 
strated experimentally in NbN superconducting wires for sin- 
gle photon detection'^ On the other hand, the geometric 
inductance of the loop Lc = 47tAl/{c^Ii) is negligible in com- 
parison with the kinetic one when A^/A » 1, where Ai is the 
area and is the length of the loop. 

The SCNW shunted by an external resistor can be repre- 
sented by an effective circuit as shown in Fig. [T] We use the 
resistively capacitively shunted Josephson junction (RCSJ) to 
model the PSCH A slightly different setup with the induc- 
tor in the external shunt branch was discussed in Ref. "76" for 
studying the JJs with an external shunt circuit and the effect 
of a parasitic inductance. In the following we provide a fully 
quantum mechanical description of the circuit in Fig. [TJb) and 
derive the equations of motion for the superconducting phase 
and the current through the PSC'^ Here, is a phase differ- 
ence across the PSC, Q is the canonically conjugate charge at 
the PSC with the commutation relation, [0, Q] = lie where e 
is the charge of an electron. The phase differences across the 
shunt resistance and the inductance are 0r and (^r -0) respec- 
tively. The capacitance of the PSC is Cpsc- The applied bias 
current is /b, and /i is an operator for the current through the 
inductance (and the PSC). It is given by /j - S(0R-0)/(2eLK). 
The full Hamiltonian reads H - Ho + Hbp + Hbs, where 



Ha = 



2C 



■cos(/> + ^— j — ^1^^ —Jb^r, (1) 
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where Hq is the Hamiltonian of the PSC plus the kinetic in- 
ductance and the coupling with the current bias. Here 1^ is 
the critical cuiTent through the PSC. We model the dissipa- 
tion arising from the quasipaiticle channel in the PSC and the 
external shunt resistor by incorporating two bosonic baths of 
harmonic oscillators. The bath related to the junction is H^p 
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and that with the shunt resistor is H^s 

d2 ,.,2 



oo pZ . 



OO pZ. ^ oo 



(2) 



The operators Xj, xj and Pj, pj are conjugate positions and 
momenta of the bath oscillators with frequencies ajj, cbj. The 
coupling between the oscillators position and the phase opera- 
tors are chosen to have an adequate description of noises from 
the baths. We choose ohmic dissipation from the baths by fix- 
ing the coupling parameters cj, Aj and frequencies atj, wj^. 
Following the procedure in Ref. [18] we employ the Heisen- 
berg equation to write equations of motion for ^, Q, Xj, xj, Pj 
and pj. We integrate out the baths' degrees of freedom from 
these equations by the exact solutions of the baths' degrees 
of freedom along with the ohmic responses from the baths. 
Thus we obtain the following generalized quantum Langevin 
equations for 4> and Ii , 

Ti dd> h (fid) 

^■^^^^S-(^B-/:)/^. + y.. (4) 
at 2e at ^ ' 

Here, and Y^^ are Gaussian current and voltage noises from 
the PSC and the external resistor The noise properties are, 
(/np) = 0, (Vns) = and the spectral function 



5„pM = fl'f— <{/„p(f)/np(0))) = — coth — — , 
J 2 Rv<^c ^IknT' 



Here is the Boltzmann constant and T is the temperature 
of the wire. The self-heating in the nanowire for a small ex- 
ternal shunt resistor is weak, and T is equal to the electrode 
temperature. Though the solution of the above generalized 
quantum Langevin equations is important and necessary for 
studying the macroscopic quantum tunneling, here we are in- 
terested in the classical dynamics of the system. Thus we re- 
place the operators by their expectation values, for example, 
^ by (p. We also choose high temperature limit of the noise 
correlations, and replace them by S„p{a>) - Ik^T/Rpsc and 
S „sia>) = Ik^TRs- The classical noise correlations in the time 
domain are given by 

</np(?l)/np(?2)> = 2kBTRplcS(tl - tl), (5) 
</n,s(fl)/„s(f2)) = 2kBTR-'6(t, - f2), (6) 

where we define y„s = RJtvi- Now we introduce the following 
dimensionless variables: /b - I^lh, h - hlh, 'nP - hvllc, 
'ns = hslh are the normalized currents, - R^/Rpsc is the 
normalized external shunt resistance, f' = cj^t is the normal- 
ized time with Wp - y/lelJihCpscj, Ik = L^cop/Rp^c is the 
normalized inductance, and T = k^cOpT / (J3RpscIc^ is the nor- 
malized temperature. Next we rewrite the classical version 



of the equations of motion Eqs. ([3| and Q and the classical 
noise correlations in Eqs.(|5]l and ^ using the dimensionless 
variables. 



dip d^(p 
ii = sm^+jS— + —r +i„p, 
dt' dt'^ 

di\ d(p 

^K-TT +j6t7 = *^'B - ll)rs + insr., 

dt dt 
(inp(fi)inp(t'2)) = 2/3T6(t[ - t'2), 
(in.(t\)ins(t2)) = 2/3r;'T6(t\-t'2), 



(7) 

(8) 

(9) 
(10) 



where /3 - ^tij {2eI^R^^^Cp^c) is the damping factor. For 
typical parameters of MoSi superconducting nanowires"^, we 
have Rps,c ~ 100 Q, Cpsc = 8 pF, = 5.55 ;uA and = 
50 GHz, jS = 0.1. As jS <«; 1, thus the dynamics of the PSC is 
underdamped. For I - 90 nm, A - 100 nm^, A = 700 nm^S 
we estimate the inductance Zk ~ 1- Therefore, the kinetic 
inductance is not negligible in the practical nanowires. We 
can find ii from Eq.(j8]l, 

iM)^h{t', =0)e-''./^ 



Jo ^ 



+ ill 



n dt 



y)dt'2. 



(11) 



with T = Ik/ fa- After substituting Eq. ( [TT| l in Eq. (|7]l and 
using /[(f'j = 0) = /b, we obtain 



dt'/ 



where 



d(p 

y(fj - t2)—dt2 + sin0 - /b + in - 0, 



(12) 



6(t[ - t'2) + 



in - inP - I 

Jo T 



-ins(t'2)dt'2. 



(13) 



The presence of the inductor introduces a memory effect in 
the damping coefficient even though we choose ohmic dissi- 
pation from the baths at finite temperature. The fluctuation- 
dissipation theorem is still valid for the nonlocal dissipation 
and the noise, (^in (f'J in (4)) - '2Ty(t[ - t'^). 

III. V-l CHARACTERISTICS AT T = 

Here we investigate the effect of the kinetic inductance on 
the V-I characteristics in the absence of noises. Equation ^ 
describes the dynamics of a phase particle in a tilted wash- 
board potential, U((p) - 1 - cos0 - When /b is ramped 
up adiabatically, the voltage induced by this slowly increas- 
ing current is negligible, and thus ii - ip,. The phase par- 
ticle always stays at the energy minimum of the potential at 
= sin"' /b until such a minimum disappears at ip > 1. 
Then the system switches to the resistive state. Therefore, 
the switching current I^^ is equal to 1^ for a slowly varying 
bias current, and I^^ is independent of Ik- The inductance Zk 
is only effective in the dynamic state, and the switching to 
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FIG. 2. (color online) V-I characteristics at T = obtained by nu- 
merically solving Eqs. jTJ and ^ without the noises for two different 
/k. Weuse/3 = 0.1. 



the resistive state is unaffected by the presence of the kinetic 
inductor at T = 0. 

The V-I curve becomes hysteretic for underdamped JJs, 
and the junction is retrapped to the superconducting state at 
/r, where the resistive state becomes unstable. The presence 
of the inductor changes the dynamics of the superconducting 
phase. For a small inductance t - /k/'"s ^ 1, we obtain from 
Eq. ([8| to a leading order in t, 

'"s at' dt"^ 

Substituting ii into Eq. (|7]i, and rescaling the time f' = 
t" -^(1 -ySr/rs), we aiTive at 



iB = sin + 



(J3+/3/r,) d(f> d^ 



dr 



^j\ -/3T/r,dt" 



,2 " 



(14) 



The retrapping cuiTent increases with the damping coefficient 
in the RCSJ model. The presence of the inductor increases 
the damping, thus the retrapping cuiTent is increased with an 
increasing inductance for a small /k <k 1 . 

We solve Eqs. (|7]i and ([8]l numerically for an arbitrary value 
of Zk- The scaled voltage v = V/(^psc^c) is calculated from 
dip/dt'. The typical V-I curves are shown in Fig. |2]in the 
absence and presence of an inductance for different shunt re- 
sistances. The zero-temperature retrapping current ;V depends 
on the inductance while the zero-temperature switching cur- 
rent isH. is independent of Ik- At a given l^, i,- increases with 
a decreasing and the hysteresis finally vanishes for a small 
enough r,. The dependence of the zero-temperature iy on Zk 
is shown in Fig. [3] For a small r,, iV = /.nv - 1 and the 
retrapping cuiTent is independent of Ik- For a large r,, the 
retrapping current first increases with an increasing Zk for a 
small Zk which is consistent with the analysis for a small Zk in 



FIG. 3. (color online) Dependence of the zero-temperature retrap- 
ping current on inductance Zk for different shunt resistances r^. 



It decreases with an increasing Zk. The presence of an in- 
ductor does not affect the zero-temperature /jw, and the resis- 
tance of the resistive nanowire. The dynamical resistance is 
= dV/dlE ~ Rpsclfiil + l/'"s)]"' in the resistive state. 



IV. STOCHASTIC SWITCHING AT T > 

At a finite temperature, upon ramping up the current, the 
system may switch to the resistive state at a bias cuiTent below 
the mean-field switching current I^v, - 4 due to thermal fluc- 
tuations. The distribution of Isw is determined by thermally 
activated escape rate of the phase particle from the meta-stable 
potential. The presence of an inductor induces a memory ef- 
fect into the dynamics. In this section, we study the effect of 
the inductance on the rate, and then calculate the distribution 
of the switching current from the rate, which enables a direct 
comparison between the theory and the experiment. 

When Zk is small, r - 1^1 <k 1, Eq. 
order in t becomes 



12 1 up to a linear 



Zb = sin0 -I- y8 + 



P_ 
Ts / df 



+ 1 - -T 



/ dr 



V/2 



(15) 



(!n(?'i)Z„(4)> = 2f L8+^|5(f'l-4). 



Equation (15 1 describes the motion of a particle with mass 
M - (\ - prlr^) in the potential U{(p), and the thermal escape 
rate F for this problem is well knowrP^, 



exp[-At//f], (16) 



Eq. ( 14 1. For a comparatively larger Zk, the zero-temperature at the barrier (well). 



w here ^ {J5 + /3/r,) AU ^ U(4>h) - U((f>o), oj,, 
^-djU(4>h), and ljq = ^d^Uicpo)- Here 4>h (0o) is the pha 
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For a general value of t, the switching rate with a nonlocal 
damping coefficient has been calculated in Ref. 21, and we 
here use their results. 



2na>b 



exp[-AU/T], 



(17) 



where a is a pole at the right-half plane of the function f(z) - 
[z2 -col+ zL[y(n]]-\ where L[y(f')] = + ZT)(r,)]-i 

is Laplace transform of yit'). The kinetic inductance does 
not change the barrier height AU, but modifies the pre- 
exponential factor of the rate. The rate is still dominated by 
the barrier, but the enhancement of the rate can be large by 
taking the inductor into account in some parameter space. For 



a small r, Eq.( 17 i reduces to Eq. (16 1, while for r — > oo 



r = 



2a)h 



In 



exp[-At//r]. 



(18) 



In the last limit, the effect of kinetic inductor and shunt resis- 
tor vanishes and the rate is the same as that in a single JJ. The 
dependence of the switching rate on /k at two different is 
shown in Fig. |4](a) and (b). The rate increases with /k. For 
/k = 0, the effect of the shunt resistor is to renormalize the 
damping of the nanowire fi^s - P + P/fs, which reduces the 
rate. When /k is quite large, the current in the wire does not 
change much with time, ii ~ and the shunt circuit becomes 
less important. Thus the switching rate reaches the limit of a 
JJ without the shunt resistance. For Rf^/Rpsc ~ 1, the incre- 
ment of the rate is small, while for a small shunt resistance 
Rs/Rp5C ^ 1, the enhancement of the switching rate can be 
large due to the presence of the inductor The dependence of 
the switching rate on r, at two different Ik is shown in Fig. |4] 
(c) and (d). The rate increases with r,, but the increment is 
small for a large Ik- 

We proceed to calculate the distribution of the switching 
current from the rate. In experiments, one sweeps the bias cur- 
rent slowly to avoid the nonequilibrium heating effect (much 
slower than any other characteristic time scale in the system), 
and measures the voltage. Once there is a nonzero voltage 
signal, a phase slip event is counted. One then resets the 
system to the superconducting state by reducing the current, 
and repeats the measurements, from which the distribution of 
switching current is constructed. The probability for the sys- 
tem leaving the meta-stable potential after time t during the 
process of ramping up the bias current is^^ 



WMf')] = 1-exp 



Jo 



r[/B(f'i)] dt\ 



(19) 



When is is increased at a constant rate ib = disldt' — 
constant, the distribution function is then given by 



PUb) 



dW 
diB 



is 



exp 



Jo 



di'fi 



iB 



(20) 



The distributions of /s„ for several temperatures are shown in 
Fig. |5] The linewidth of the distribution increases with an 
increasing temperature. The distributions in the presence and 
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FIG. 4. (color online) (a) and (b): Normalized switching rate 
T{Ik)IT{Ik = 0) as a function of inductance Ik for two different shunt 
resistance r,. (c) and (d): Normalized switching rate r(rj)/r(r, = oo) 
as a function of for two different Ik- The normalized switching 
rates are temperature independent. 



absence of an inductor are shown in the inset of Fig. [5] The 
shape of the switching current distribution is modified in the 
presence of the inductor. 




FIG. 5. (color online) Distribution of the switching current at 
different temperatures f . The inset shows the distributions of 
with and without the inductor. The parameters for the main panel are 
r, = 0.5,/? = 0.1, /jf = 4 and ig = 0.01, and for the inset is r, = 0.02, 
P = 0.2, t = 0.04 and Ib = 0.01. 
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V. DISCUSSION 



In short, we have proposed a new effective circuit for ex- 
ternally shunted superconducting nanowires by taking the ki- 
netic inductor of the superconductor into account. The pres- 
ence of inductor introduces a nonlocal damping in the dynam- 
ics of superconductivity, (a) First it minimizes the heating 
in the nanowire by introducing a delay in the redistribution 
of current between the resistive shunt and the nanowire. (b) 
The kinetic inductor does not affect the mean-field switching 
current, but it modifies the mean-field retrapping current in 
a non-trivial way. (c) We also consider the thermally assisted 
switching from the superconducting state to the resistive state. 
The switching rate increases with the inductance. The result- 
ing distribution of the switching current is also modified. 

Inductive shunting plays a significant role in building su- 
perconducting circuits using Josephson junctions^"*. One of 
the interesting direction for potential application of SCNWs 
is to use them as superconducting qubits in quantum infor- 
mation processing and circuit quantum electrodynamics type 
set-ups. Therefore, the inherent kinetic inductance of SCNWs 
may provide an advantage for using the nanowires compared 



to the Josephson junctions. A voltage biased narrow SCNW 
with coherent quantum phase slips is exact dual to a current 
biased RCSJ modePl It has been argued recently that weak 
links with slightly larger resistivity occur naturally in appar- 
ently homogeneous nanowires, and they can localize quantum 
phase slips by quantum fluctuations at low temperatureJ^^l 
Here we show that the phase slip center in series with an 
inductor is responsible for the dynamics in the externally 
shunted SCNWs at temperatures close to the critical. Since 
the inductance depends on the length and the cross section of 
the nanowire, the role of the kinetic inductor can be tested 
in experiments by measuring the voltage-current characteris- 
tics and the distributions of the switching current in resistively 
shunted nanowires of different sizes. 
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